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DYNAMIC MODEL OF TEE EAHIB'S UPPER ATMOSPHERE 


1.0 INTRODUCTION 

Although high resolution data from satellite— borne mass spectrometers, 
accelerometers and other instruments have contributed greatly to our under- 
standing of the phy'sics of the neutral upper atmosphere, the older atmos- 
pheric models, based entirely on satellite drag analysis, remain viable in 
most applications to orbital mechanics. The more recent models, such as 
MSIS and 311 , represent composition and changes in composition more accu- 
rately than do the models based on drag. For the most part, however, they 
are not significantly better in their representation of the total density. 
This, combined with the greater complexity and computing costs of the more 
recent models, argues strongly for continued use of the older models in 
those situations where the representation of total density and variations 
in total density are the only criteria of model performance . 

The drag— based models are deficient in the way they represent the 
large variations in total density that are associated with geomagnetic dis- 
turbances, however. These variations are represented in the models as 
being uniform over the globe, while the more recent high resolution data 
have shown them to have very significant global structure, mainly in rela- 
tion to geomagnetic latitude. This is a matter of concern in a number of 
currently important applications where it is necessary to accurately 
predict the density and its variations on a local basis. For this reason, 
we undertook the task of modifying the J70 atmospheric model used at MSEC 
to include a high resolution model of the geomagnetic variation in the 
thermosphere and exosphere. 

The model of the geomagnetic variation that we included is an improved 
version of the one that we developed from analysis of the ESRO 4 mass 
spectrometer data and that was incorporated in Jacchia's 1977 comprehensive 
model. It is described in detail in Appendix A. It includes the variation 
with geomagnetic local time as well as with geomagnetic latitude and also 
includes the effects due to disturbance of the temperature profiles in the 
region of energy deposition. Tie implementation of the model in the MSFC 
program is described in the next section. 
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2.0 MODEL Iitf>LEMENTATION 


The high resolution model of the geomagnetic variation desoxited in 
Appendix A was implemented in computer software by modification of the 
J70—3 program used at MSEC. A copy of the FORTRAN source file for this 
program was supplied to us in card form. The source file for the modified 
program, which we designated J70-3X, has, in turn, been shipped to MSEC in 
card form together with a compiler listing of the program and samples of 
the output that it produces. 

Modification of the program included revision of the main program and 
the TINF and JAC $Dbroutines and addition of the subroutines GANG, CONGI 
and AVER. The labeled common block MAGMOD was introduced to aid in commu- 
nication between the different program units. The hyperbolic tangent was 
also introduced in several locations as an assumed library function. 

The modified program provides for optional calculation of the geomag- 
netic variation according to either the 1970 Jacchia model, as in the 
original program, or the 1971 Jacchia model rather than the high resolution 
model. In addition, the high resolution model may be specified in either 
of two diurnally— averaged forms rather than in unaveraged form. 

21 Modes of Operation 

The modified program operates in one of five possible modes relative 
to the calculation of the effects of the geomagnetic variation. The par- 
ticular mode is specified by a parameter that was added at the end of the 
"case" list that is a part of the control input to the program. The 
addition of this parameter, MQKi, was the only change made to the program 
input. Calculation of the geomagnetic variation in the different modes is 
described briefly in the following subsections . 

2.1.1 571 Geomagnetic Model. MCI® » —I 

In this mode, the geomagnetic variation is calculated according to 
equations 18 and 20 of Jacchia's 1971 model atmosphere . Equation 18 
provides only for an increase in exospheric temperature and applies to 
heights above 350 km. It is identical to the equation given for the geo- 
magnetic variation in Jacchia 's 1970 model and was derived in one of our 
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earlier studies of satellite drag. It does not give good results below 
about 200 Ion because the variation of density with exospheric temperature 
in the static models is too small to adequately represent the geomagnetic 
variation at those heights. Equation 20 includes a direct increase in 
density in addition to an increase in exospheric temperature and gives 
better results at low heights, We employed the same interpolation scheme 
used in the implementation of the 1971 model included in the 1972 CIRA to 
make a smooth transition between the two equations in the vicinity of 350 
km. This mode was included in order to provide a more suitable standard of 
comparison for the high resolution model at low heights, 

2.1.2 J70 Geomagnetic Model. MIE = 0 

In this mode, the geomagnetic variation is calculated according to the 
equations given with Jacchia's 1970 model atmosphere. Results in this mode 
are identical to those obtained with the original version of the program. 
Note that the program provides for calculation using either equation 21 or 
equation 22 of the 1970 model, depending on whether the Kp geomagnetic 
index or the equivalent &p index is specified in the input. At present, 
however, the relevant indicator is set internally in the program so that 
the value in the input is always assumed to be the ap index.* 

2.1.3 Local —Time —Averaged Model . MdE = I 

In this mode, the high resolution model of the geomagnetic variation 
is used, but an average is taken over local time. The results represent 
daily averages for points above a fixed location on the surface of the 
earth and will depend on the geomagnetic latitude that corresponds to the 
geocentric latitude and longitude specified in the input . This mode was 


* This is true regardless of the operating mode that is specified. None 
of the other modes provide for direct calculation using the index, how- 
ever. Instead, the Sp index, when indicated, is converted to the equiva- 
lent Kp index by a call to the subroutine CONGI prior to entry in the 
programmed equations. 

It would have been appropriate to have the program print the geomag— 
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provided mainly for diagnostic and developmental purposes. 

2.1.4 Longitpdc~Avor»ged Model s MIX » 2 

In this mode, the high resolution model of the geomagnetic variation 
is used, but an average is taken over longitude. The apparent solar time, 
as determined by the geocentric longitude and the GMT given in the input, 
is kept constant. In this case, the results represent daily averages as 
would be seen by an orbiting object at the point in its orbit corresponding 
to the geocentric latitude and longitude and the GMT specified in the in- 
put. This mode would, therefore, be appropriate in orbital integration or 
a comparison with drag measurements where the time step is a day or more. 
Note that, while considerable smoothing will occur in high latitudes, a 
strong local time dependence will remain in middle and low latitudes i n 
addition to the dependence on latitude. For this reason, the program 
prints the apparent solar time in the output along with the specified lati- 
tude and longitude. Also note that both diurnally— averaged modes (MODE = I 
or 2) are somewhat more expensive to run than the non— averaged mode (MODE = 
3), since the averages are obtained numerically* by summing values taken at 
10— degree intervals around the globe. Although .no detailed comparisons 
were made, this is not expected to be very significant compared to the 
total cost of running the program, however. 

2.1.5 Complete High-Resolution Model, MOIX » 3 

In this mode, the geomagnetic variation is Calculated according to the 
improved high— resolution model developed from mass spectrometer observa- 
tions. Most of the details of this model, including the values of the 
constants used, are given in Appendix A. In this mode there is no averag- 
ing and the results are directly dependent on both the geomagnetic latitude 
and the geomagnetic local time, both of which are calculated from the nor- 


netic latitude in the output along with the specified geocentric latitude 
and longitude in this mode. This, however, was overlooked in the program- 
ming process. 

• The averaging is done in subroutine AVER for either diurnally— averaged 
mode m 


i 


4 



mal input to the program. The calculation is that for a centered dipole.. 
The geomagnetic latitude and geomagnetic local time are both printed with 
the specified geocentric latitude and longitude in the output. As in all 
of the modes involving the high resolution model, the temperature profiles 
are assumed to be perturbed according to equation 32 of Jacchia's 1977 
atmospheric model, 

AqT(z) = AgT® tanli[0.006 (z - 90.0)] , 

where z is in km. This equation is incorporated directly in the integra- 
tion of the diffusion equation in subroutine JAC. It permits the effects 
of geomagnetic disturbance to be accurately represented even at very low 
thermospheric heights. 

2.2 Comsotatioa of Geomagnetic Coordinates 

The high resolution model is a function of both the geomagnetic lati- 
tude and the geomagnetic local time. These angles are computed from the 
geocentric latitude and longitude and the solar declination and hour angle 
in subroutine GANG, assuming the earth's magnetic field to be represented 
by a centered dipole. 

The geomagnetic latitude, 0G» is easily compated''' from 

sin0g= sin 0 sin 0p + COS0 cos 0p cos(X - Xp) , (1) 

where 0 and X are the geocentric latitude and east longitude, respectively, 
of the point in question and 0p and Xp are the geocentric latitude and east 
longitude, respectively, of the boreal geomagnetic pole. We assumed 0p = 


* The (adiabatically) invariant geomagnetic coordinates seem to correlate 
best with small details of the geomagnetic variation, but the differences 
between the two coordinate systems are small and the dipole coordinates are 
entirely adequate for our purposes here. 

The coordinate conversion described in this section is illustrated in 
"Coordinate systems based on the geomagnetic field " in the International 
Dictionary of Geoohvsics . Vol. 1, pp. 616-617. The equations given in that 
article are not fully developed for machine computation as are those given 
here, however. 
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79.0 degrees and Xp = 291.0 degrees. 

To compute the geomagnetic local time, we first compute the "|eooag- 
netic longitude", Xg, which is the angle, measured counterclockwise at the 
geomagnetic pole, between the geographical meridian extending southward 
from the geomagnetic pole and the great circle through the geomagnetic pole 
and the point in question. We have 


cos 0 sin(X-X|j) 
s i n k = -r 

COSVQ 

sin 0 - sin 0(5 - sin 4 

cosXq = 

cos 0p cos 0Q 


( 2 ) 


These equations become Undefined in the case where cos0q = 0 (i.e., where 
the point in question is at one of the geomagnetic poles). In that case, 
we set Xq = 0. 

Next, we compute cos p and sinp, where p is the angle between the 
geomagnetic pole and the direction of the anti— sun. We have 


cosp = “ s in0p s in 8 g + cos0p cosSg cos( t +Xp ” X) 
sinp = ( 1 -cos* P)i', 


(3) 


where 85 is the solar declination angle and t is the apparent solar time {= 
solar hour angle + n radians) at the point in question. Taking sinp to be 
always positive merely implies that P lies between 0 and n radians, as it 
should. 

We then compute the "geomagnetic time " of the geographical meridian 
that passes through the geomagnetic pole. This we define to be the angle, 
measured counterclockwise at the pole, between the great circle through the 
pole and the anti-sun position and the southward extension of the meridian 
that passes through the pole. Calling this angle Atg, we have 


s in Atg = 


cos 8 g sin(t + Xp- X) 
sin P 


cos Atg = 


sin 85 + cosp sin 0 p 
siiiP cos 0 p 


(4) 


Since neither sinp or cos 0p can be equal to zero, there is no problem with 
these equations becoming undefined. We can then compute the geomagnetic 
local time, tg» which is given by 
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to = X,6 + At(j. 

A good approximation to the geomagnetic local time is given by 


(5) 


tG»t + Xp + XG-X, (6) 

which eliminates evaluation of equation 3 and the determination of Atg from 
equation 4. The savings that could be gained by using the approximation 
were not thought to be significant in terms of the total cost of running 
the program, however, and exact computation of the geomagnetic local time, 
utilizing equations 2 through 5, was implemented. 

2.3 Eauatorial Wave Modification 

As given by equation A-12 of Appendix A, the "equatorial wave " compo- 
nent of the high resolution model consists of a proportional increase in 
density that is independent of height. That this component of the geomag- 
netic variation should be independent of height in the lower thermosphere 
is not physically realistic and represents a flaw in the formulation of the 
model. To correct this flaw in the implementation of the model, we devel- 
oped the following expression to replace the original equation for the 
equatorial wave: 

Ae log nj =Ae log p = 5.55 1 10“ * tanh[0.006(z - 100.0)]Acos» 6^ (7) 

where p is the total density and A is the polar increase in exospheric 
temperature given by equation A-3 of Appendix A. With this modification, 
the high resolution model gives density variations at low heights in the 
equatorial region that are in good agreement with those reported by Marcos 
£± al * from lor— altitude accelerometer measurements made in a number of 
satellite experiments at heights down to about 140 km. 


* Marcos, FA., DJF. Gillette, and E.C. Robinson, 1982: A Global Thermo- 
ssheric Density Model Based on Satellite Accelerometer Data. AFGL-TR-82- 
0025. 
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24 Weight^.d— Mean Geomaxs 0 tio Index 


The geomagnetic index is specified as one of the input parameters to 
the J70-3 atmospheric model program. TiiBS, we were not able to implement 
calculation of a weighted— mean value of the index as specified by equation 
A— 5 of Appendix A If the program is to be used to simulate actual varia- 
tions in the atmosphere, we strongly recommend that this be done, however, 
and that the weighted— mean-index be used as input to the program.* We would 
be happy to assist MSPC in whatever way is required to achieve this end. 

Since we first introduced the idea of using a weighted— mean index in 
our model of the geomagnetic variation, we have discovered confirmation of 
it in the work of others in the field. Nisbet e t al '*' combined this idea 
with use of an auroral electro jet index in a study of atomic oxygen varia- 
tions observed by the OSS mass spectrometer on the AE-C satellite during 
two geomagnetic storm periods. They obtained significantly better 
representations of the observed variations than could be obtained using 
individual values of the 3— hour index. Of course, this may have been 
due in part to the superiority of the electro jet index as an indicator of 
disturbance in the atmosphere. It is worth noting, however, that the time 
constants («24 hours) that were found to give the best results for the 
recovery phases of the storms are entirely consistent with the one that we 
arrived at quite independently from our examination of the ESRO 4 and AE-C 
mass spectrometer data. More recently, Hedin^ has included an index based 
on a similar "summation formula" for use with the model of the geomagnetic 
variation in his latest version of the MSIS thermospheric model. He gives 
an example in which this index is utilized that shows how the "rise and 


* The concept of a weighted— mean of the geomagnetic index as an index of 
atmospheric disturbance is not, of course, limited to the high resolution 
model. It applies to any model of the geomagnetic variation. 

Nisbet, J.S., C, Stehle, and E. Bleuler, 1983: Initial tests of an 
index based on AL values for modeling magnetic storm related perturbations 
of the thermosphere, J. Geoohvs. Res. . 88 . pp. 2175-2180, 

^ Hedin, A,E., 1983: A revised thermospheric model based on mass 

spectrometer and incoherent scatter data: MSIS-83, J. Geonhys. Res. . 88 . 
pp. 10170-10188. 
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fall of the density predicted using the summed three hour indices is mtich 
more representative of the data than either the prediction based on the 
daily index or the variation of the 3— hour index itself." He uses separate 
time constants for the exospheric temperature and each of the various con- 
stituents and finds the values of these to vary between 6.9 hours and 13.0 
hours. His value for atomic oxygen, rhioh is most important as far as 
total density is concerned at most heights, is 11.4 hours. Ties, while he 
confirms the applicability of the concept, his estimate of the time con- 
stant is somewhat smaller than was indicated by our work and that of Nisbet 
et al . We hope to be able to pursue the question of the best value of the 
time constant further, utilizing a larger portion of the mass spectrometer 
data that we have. 
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APPENDIX A 


A NEW MODEL OF THE GEOMAGNETIC VARIATION 
IN THE IMTO ATMOSPHERE* 


A. 1 INTRODUCTION 


Large variations in the density of the thermosphere and exosphere in 
association with geomagnetic disturbances were first detected by Jacchia 
(1959, 1961) in the drag on artificial satellites. Densities obtained from 
the analysis of drag were then used extensively to study and model this 
geomagnetic variation in the upper atmosphere. Model development based on 
drag analysis culminated with the model derived by Jacchia et (1967) 
that was subsequently incorporated, with only slight modification, in the 
1972 CIRA reference atmosphere and other general models of the upper atmos- 
phere derived from drag analysis. 

The models of the geomagnetic variation based on drag analysis were 
extremely simple in form and usually represented the variation as being 
uniform over the globe. As it turns out, the geomagnetic variation is 
extremely complex in form but, while this was certainly suspected, the drag 
data revealed few details because of their poor resolution. It was not 
until data from situ accelerometer and mass spectrometer experiments 
became available that the geomagnetic variation could be modeled in any 
detail. The process is still far from being complete, however. 

There are, at present, two models based on high— resolution data in 
use. One of these is that of Jacchia et (^976, 1977) that is incorpo- 

rated in Jacchia' s most recent general model of the thermosphere and exo- 
sphere (1977). The other is that of Hedin ej al that is incorporated in 
the MSIS model (1977a, 1977b, 1979). Here, we will describe initial ef- 
forts directed toward extending and improving the model of Jacchia et al . 
In this, we have four main aims. These are: 

1) To incorporate the variation with local magnetic time. As it is 
now, the model includes only the mean variation with latitude. 
The local time variation is extremely complex, but the main fea- 
tures are quite significant and can be included without too much 
difficulty. 


• This is a reproduction of a paper originally presented at the Ninth 
Conference on Aerospace and Aeronautical Meteorology in Omaha in June, 
1983, and later published as USAF Report AFGL-TR-83-0253 , 
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To allow the form of the variation to vary with the level of 
disturbance. Slowey (1981) has shown the latitudinal form of the 
variation to broaden in high latitudes with increasing levels of 
disturbance. This effect appears to be related to a shift in the 
location of some of the heat input. Yeng (1982) has reported that 
the polar cusp region shifts significantly toward lower latitudes 
for larger disturbances, 

3) To allow for the prior heat input in establishing an index for the 
level of disturbance. The Kp geomagnetic index admittedly leaves 
much to be desired as an indicator of disturbance in the atmos- 
phere. It is increasingly clear, however, that there is much to 
be gained by taking account of the persistence of the effects of 
disturbance. Theoretical calculations such as those by Faller— 
Rowell and Rees (1981) and the semiempirical models of Hedin et al 
(1981) indicate that the disturbance history over the preceding 5 
days is probably significant . 

4) To correct what is now seen as a rather serious flaw in the 
original model in not permitting an increase in exospheric temper- 
ature at the equator. The effect at the equator was modeled en- 
tirely in terms of an "equatorial wave" of density that was as- 
sumed to originate in high latitudes and progress towards the 
equator. While it is clear from the mass spectrometer data that 
such a wave is a significant component of the low latitude 
response, it is now equally clear that a thermal component also 
exists in low latitudes. Separating the two components in data 
from a single height or from a limited range of height, which was 
all that was available to us until recently, is extremely diffi- 
cult, however, and a final determination of the relative impor- 
tance and extent of the equatorial wave will depend on future 
analysis of data from lower heights. 


A2 NATURE OF TEE VARIATION 


In figure A—1 we show a plot of molecular nitrogen density measured 
during a portion of one orbit of the Atmosphere Explorer-C satellite by the 
OSS mass-spectrometer . The data were collected during a period of rela- 
tively high geomagnetic activity (1.5 hours earlier, was 6+) and when 
the satellite was in a nearly circular orbit. Since the height varied so 
little, in was safe to reduce them to a single height for purposes of 
analysis and that is the way they are shown here. The abscissa is the time 
in minutes, starting from an arbitrary zero point. The scale of the or- 
dinate of the Nj number density is not shown, but it represents the change 
in the logarithm of the number density. The maximum density here corre- 
sponds to a change of 0.5 in lag,, n(N,), an increase by a factor of more 
than 3 in number density, relative to the density at the origin of the 
plot. The smooth curve at the top of the figure is the (adiabatic invari- 
ant) geomagnetic latitude plotted on the scale at the left of the figure. 
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Figure A—1. Logarithm of N, density at a height of 250 km as 
measured ' by the OSS mass spectrometer on AE—C during a portion of 
orbit 4919. The abscissa is the time in minutes from an arbi- 
trary zero point and the ordinate is not given (see text). The 
smooth curve is the invariant geomagnetic latitude with ordinate 
scale on the left. 

The strong dependence of the geomagnetic variation on geomagnetic lat- 
itude is well established and we can see it in orbit after orbit of the 
AE-C and other satellite— borne mass Spectrometer data, not only in If but 
in all constituents . Not all constituents behave the nay N, does, however, 
and this is another important aspect of the geomagnetic variation. The 
behavior of the lighter constituents, such as He, is just the opposite of 
that of Nj : their number densities decrease where those of Nj increase. 

And the density of atomic oxygen, which is of primary importance as far as 
total density is concerned throughout most of the thermosphere, may either 
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increase, decrease or remain unchanged. In addition to a purely thermal 
effect, there is another factor involved in which the molecular mass is 
very important. We shall represent this factor in terms of a variation in 
the height of the homopause. 

The variation is also clearly not just dependent on the geomagnetic 
latitude. There are two distinct maxima in the data of figure A-1 that 
represent separate regions of enhancement in the vicinity of the auroral 
oval. As the figure demonstrates, there is also an appreciable over— all 
asymmetry with respect to geomagnetic latitude. 


A. 3 FORM OF THE MODEL 


The model of Jacchia ^ represents the density and composition 
changes that occur in association with a geomagnetic disturbance by a local 
increase in exospheric temperature and a proportional increase in the 
height of the homopause. The increase in the height of the homopause is a 
convenient device by which to represent an effect (the inverse variation of 
the lighter constituents) that may actually be due more to wind— induced 
vertical diffusion. Superimposed on these two effects is the "equatorial 
wave" referred to in the introduction, in which the number densities of all 
constituents increase in the same proportion and that is centered on the 
equator. We will adopt the same representation here and model the change 
in the logarithm of the number density of the species i as the sum of three 
components : 


A(3 Iogn£ = Aj lognj + Ag logn^ + Ag lognj (A-1) 

where Aj log n^ is the thermal component, Ag log nj is the component due to 
the change in the height of the homopause and A ^09 Oi -is the component 
due to the equatorial wave. The thermal component Aj log nj is to be 
evaluated from an atmospheric model assuming an increase in exospheric tem- 
perature AqTcd given by 


AgT» = AF(0,X) (A-2) 

where A is the amplitude given by 

A =57.5K^[l + 0.027 exp(0.4Kp] , (A-3) 

as in the earlier model, and F(0,X) is given by 
F(0, X) = Box + “os sin“0 + 

cos* 0(aji sinX + axa cosX) + 
sin20(axx sin X + a* a cos X) + 

sin'* 0 siii20(a, ^sinX +aaa cosX) , (A-4) 
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where t and X. are the geomagnetic latitude and local magnetic time, respec- 
tively. To take the effects of persistence into account, we assume Kp' in 
equation A— 3 to be the weighted mean of the lagged 3— hourly Kp geomagnetic 
index taken over the 41 values in the preceding 5— day interval as follows: 


^U) 


t£=t— c-5 

t-T , . 

ti=t-x-5 


(A-5) 


where t is the time in days, c = 1.0 d~^ and t is the time lag given by 

X =0.05 +0,1 cos* 0 (day) (A— 6) 

In equation A— 4, the first two terms represent the mean latitudinal varia- 
tion of the increase in exospheric temperature. The remaining three terms 
represent the variation of the temperature increase with local magnetic 
time (LMT) in low, middle and high latitudes, respectively. To introduce a 
variation in the mean shape of the latitudinal variation with the level of 
disturbance, we take a in equation A— 4 as 

a = 5.0-Kj^/3.0 (A-D 

A preliminary determination of the remaining parameters in equation A— 4 was 
made by a least squares fit to values of the exospheric temperature in- 
crease inferred from Nj measurements made by the ESR04 mass spectrometer. 
Exospheric temperatures were obtained from densities by inverse interpo- 
lation in an atmospheric model. The temperature increase was taken to be 
the difference between the exospheric temperature corresponding to the ob- 
served Nj density and that corresponding to the density computed for quiet 


conditions from the ES804 "quiet time 
resulting coefficients are 

" model 

of von Zahn _et (1977) . 

The 

aoi = .1425E+00 

ap 2 = 

.8137E+00 

(A-8) 

an = .1184E+00 

&12 = 

-.3604E-01 


an = -.7354E-01 

822 - 

.1038E+00 


831= .3706E+00 

832 — 

-.1441E+00 



The isotherms of the relative temperature increase as given by equa- 
tion Ar^ in the case where = 3 (n = 4) are plotted over the globe in 
geomagnetic coordinates in figure A— 2. As can be seen, there is a consid- 
erable asymmetry with respect to the pole. The maximum is between 6^ and 

9h 

at a latitude of about 80 degrees. This probably reflects both Joule 
heating by the westward auroral electro jet and particle precipitation in 
the cusp region. The extent of the asymmetry can perhaps be better seen in 
figure A— 3, where the profiles for 6^ and 18“ LMT ” near the extremes ~ are 
plotted together with the mean latitudinal profile. Of course, the differ- 
ence in high latitudes as seen by an orbiting object would tend to be 
smoothed out over intervals on the order of a day because of the earth's 
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rotation. This would not be the case in low latitudes, however, and it is 
interesting to note that, even with the poor resolution of satellite drag, 
Roemer (1971) was able to detect this asymmetry in densities dedocefl from 
drag. 



Figure A— 2. Global isotherms of the relative temperature increase 
due to geomagnetic disturbance as given by equation A-4 for = 
3 m The coordinates are geomagnetic latitude and magnetic local 
time. 


We should point out that a geomagnetic disturbance will alter the 
atmospheric temperature profiles, so that the procedure we have outlined of 
entering exospheric temperatures in a static model to determine the effects 
on the number densities is not strictly valid. There will be a height 
dependence that we have not accounted for that becomes an important consid- 
eration at lower heights. In this connection, we have developed disturbed 
temperature profiles that give good results in representing observed den- 
sity variations at heights as low as 150 Tua. We have also developed ex- 
pressions to represent the effects of these disturbed profiles on the num- 
ber densities analytically. These expressions are rather involved, so we 
will not repeat them here, but they represent an integral part of the 
thermal component of the geomagnetic variation as we would model it. 

Concerning the other two components of the geomagnetic variation, we 
would include them in the form given recently by Slowey (1983), Specifi- 
cally, we would compute the component due to the change in the height of 


Ar-6 

i 




Figure A-3. Latitude profiles of the relative temperature in- 
crease for = 3 as given by equation A-4 for 6^ and 18*^ LWT, 
The profile of the mean latitudinal variation is also shown. 

the homopause from 


Ag lognj = aj Az^, 

where Azg (meters) is to be computed from 

Azg — 22 .0 AqTq, , 

where A<jToo is given by equation h-2. and the ai's are; 

a(Ar) = + 3 ,07 X 10~* (mks) 
o(0j) » +1.03 X 10“* (mks) 

a(Na)= 0.0 

a(0) *“5.75xl0~* (mks) 
o(He) = “ 6 .30 X 10-* (mks). 

The component due to the equatorial wave would be computed from 

Ag log nj - Ag log p = 5.2xlO"*A cos* 0 


(A- 9) 


(A-io; 


(A-ll) 


(A-12) 
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where p is the total density and A is given by equation A— 3. As was 

mentioned in the introduction, some work remains to be done to better 
separate this component from the thermal component in low latitudes . 
Suitable data are available from both mass spectrometer and accelerometer 
experiments . 
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